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ABSTRACT 

Dark matter sub-halos create gaps in the stellar streams orbiting in the halos of galaxies. We 
evaluate the sub-halo stream crossing integral with the guidance of simulations to find that the linear 
rate of gap creation, 72.u, in a typical Cold Dark Matter (CDM) galactic halo at 100 kpc is TZ\j ~ 

0.0066Mg"°-^^kpc"^Gyr"\ where Ms{= M/IO^M©) is the minimum mass halo that creates a visible 
gap. The relation can be recast entirely in terms of observables, as TZ^ ~ O.OSOui"*^-^^ kpc~^Gyr~^, 
for w in kpc, normalized at 100 kpc. Using published data, the density of gaps is estimated for M31's 
NW stream and the Milky Way Pal 5 stream. Orphan stream, and Eastern Banded Structure. The 
estimates of the rates of gap creation have errors of 50% or more due to uncertain dynamical ages and 
the relatively noisy stream density measurements. The gap rate-width data are in good agreement 
with the CDM predicted relation. The high density of gaps in the narrow streams require a total halo 
population of about 10'°^ sub-halos above a minimum mass of 10^ M©. 

Subject headings: dark matter; Local Group; galaxies: dwarf 



1. INTRODUCTION 

The dark matter halos in which galaxies are embedded 
are satisfactorily modeled as smooth quasi-isothermal 
spheroids, such as Hernquist (1990) or Navarro-Frenk- 
White (NFW, 1997) f unctions, as describe d in au thor- 
itative texts such a s iBmnev fc Merrifieldl ()1998f ) and 
iBinnev fc Tremaind (|2008[) . Such mass models, with 
addition of the visible stars, gas and dust, can ade- 
quately account for most of the observed internal kine- 
matics of galactic systems. On the other hand, N-body 
simulations of the formation of dark matter halos ro- 
bustly predict that approaching 10% of the mass of any 
galactic dark halo should be in the form of sub-halos 
with numbers that rise steeply toward s lower masses 
jKlypin et al.' '1999'; 'Moore et al.' '1999"; iSpringel etahl 
[2008; Dicmand. Kuhlen k Madau 2007). The baryonic 
disk and bulge components of a galaxy will lead to sub- 
halo depletion in the central regions (jD'Onghia et al.l 
|2010[) relative to a pure dark matter model. Although 
galaxies do contain visible dwarf galaxies embedded in 
dark matter sub-halos, the numbers are far short of the 
many thousands expected. Hence, it appears that either 
the predicted sub-halos are very dark, or, substantially 
not present. 

An accurate census of the relatively low mass com- 
pletely dark sub-halos remains a significant challenge. 
In principle the dark matter itself must have some cross- 
section for interaction with either itself or baryons, which 
could potentially give rise to annihilation radiation, but 
so far there is no clear association of gamma-rays with 
known sub-halos ([Bucklev fc Hoopeil 120101 : iVivier et al.l 
1201 It) . The gravitational effects of the large num- 
bers of sub-halos quickly average with distance. Sub- 
halos are potentially detectable throu gh micro-lensini 
of muhiple- image background sources (jMao fc Schneide: 
119981: IDalal fc Kochanek|[200l . The anomalous flux ra- 
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tios of multiple image gravitational lenses are explained 
with the addition of at least one moderately massive 
sub-halo projected onto the central region of the galaxy 
(jVegetti et al.ll2010l ). but this is not a detection of a vast 
population. Over a Hubble time sub-halos lead to heat- 
ing of some ~ 30kms"^ ([Carlbcrg 2009). Such heating 
is not easily detectable in most of the visible components 
of a galaxy. 

The ongoing discovery of thin, low mass, stel- 
lar streams in galaxies opens up the possibility 
that the sub-halos can be detected through their 
gravitational effects on star streams. In partic- 
ular, sub-halos that pass near or through a star 
stream create a gap in the s t ream (llbata et al.l 120021: 
Siegal-Gaskins fc Vallurll 120081: lYoon. Johnston fc Hoed 



20111: ICarlberg et al.l 120111: IHelmi et al.ll2011D . Detect- 
ing gaps in stream pseudo-images is observationally less 
expensive than measuring stream kinematics. The chal- 
lenge of the image analysis is that the stars in any im- 
age are overwhelmingly in the foreground or background 
and the signal above the statistical noise at any point 
in the stream is often fairly low. Detections of gaps 
and lumps have been claimed with good statistical con- 
fidence for the Pal 5 strea m (jOdenkirchen et ahl 120031 : 



iGrillmair fc Dionatosll2006D . the 100 kpc NW stream of 
M31 (ICarlberg et aLlboTTI). the E astern Banded Struc- 
ture of the galaxyjGrillmaii"2011') and discussed for the 
Orphan stream (Newbcrg et al. 2010). 

The purpose of this paper is to develop a basic analytic 
prediction of the rate at which sub-halos create visible 
gaps in star streams as a statistical measure of the num- 
bers of sub-halos present. The relation can be applied 
to either simulations or observational data. In principle, 
the density of dark matter sub-halos (and other bound 
structures) in a galactic halo is simply measured with 
the number of t imes they cross a stellar stream to cre- 
ate v isible gaps (|Carlberdl20"09l : IYoon. Johnston fc Hog^ 
1201 It) . The rate at which gaps are created in a stream 
is the product of the density of sub-halos, the distance 
between sub-halo and stream which creates a gap, and 
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the velocity at which sub-halos cross the stream. To use 
this simple measure of sub-halo density we need to allow 
for the mass spectrum of sub-halos and determine what 
range of velocities and impact parameters create visible 
gaps. In this paper we use a set of restricted n-body 
simulations to help develop a gap creation rate formula. 
The last step is to compare the rate predictions to the 
available data on gaps in the star streams in the Milky 
Way and M31 as a test of the CDM sub-halo predictions. 



The rate at which gaps are created in a long thin stellar 
stream is the rate at which sub-halos of mass M with lo- 
cal volume density N{r, M) cross the stream, integrated 
out to the maximum impact parameter, 6, which creates 
a detectable gap. Full n-body simulations of the sub-halo 
content of halos find that the distribution of sub-halos, 
M{r, M), is separable into a radial dependence, n(r), and 
a mass function N{M ), such that 7V(r, M) = n{r)N{M) 
(|Springel et al.ll2008| ). The N{M) function is expected 
to have a large v ariation from one galaxy to another 
(jChen et al.l 1201 ill . The rate per unit length at which 
sub-halos cross the stream to create gaps is therefore. 



7^u(r) = 




{r)N{M)v^f{v^)TTdhdvx_ dM. 



(.1) 

The factor of tt arises from the consideration that at a dis- 
tance h from the stream, for any point on a circumference 
27r6, half the stars will be headed inward in a velocity dis- 
tribution without flows. This analysis will assume that 
the velocity distribution is Gaussian and isotropic, but 
more general distributions could easily be incorporated. 

The gap rate integral. Equation [Tl depends on the 
visibility of gaps with the distance of closest approach, 
mass and velocity of a sub-halo. These factors can be 
evaluated using numerical simulations which determine 
the average outcome of a single sub-halo passing near 
or through a stream of particles. Here a star stream is 
idealized as a ring of particles on a circular orbit. The 
mass, velocity, angle of approach and impact parameter 
are all varied. After the sub-halo has passed we measure 
the height, length and density dip in the perturbed re- 
gion, from which we develop a model which can predict 
bmax in terms of sub-halo mass and encounter velocity. 
With this numerical guidance in hand, we will integrate 
over sub-halo velocities and the mass distribution, yield- 
ing the rate in terms of the minimum sub-halo mass, M, 
that creates a visible gap. A related calculation gives the 
gap lengths in terms of sub-halo masses and can be used 
to create a relation between gap densities and stream 
width, both of which are measurable in imaging observa- 
tions. 

2.1. Simulations 

An extensive set of restricted n-body simulations 
guides our approximate evaluation of the linear gap 
rate. Equation [l] in a general context. The background 
potential is a n NFW density profile as fitted to the 
iSpringel etahl (|2008D simulation of a Milky Way or M31 
like galaxy. The stream is idealized as a circular ring 
at either 30 or 60 kpc of 10^ mass-less particles. The 
ring of particles is initiated on tangential orbits at the 




Fig. 1. — The top and side views, left and right, respectively, of 
rings that a sub-halo of mass 10^ Mq, with a scale radius 0.83 kpc, 
passes near. The left column is in boxes of 51 kpc on a side showing 
the segment of the ring where the sub-halo encounter occurred. 
The right column of sub-panels are 132 kpc wide and 0.6 high, 
that is, vertically stretched a factor of 220. In all cases the sub-halo 
motion is radial. The top two rows both have impact parameters 
of 0.06 kpc, the top being motion parallel to the orbital plane and 
the middle at an angle of 30° to the plane. In the bottom row 
the closest separation is 1.5 kpc above the plane. The particles are 
shown at 6.92 Gyr, about 3.8 rotational periods after the encounter. 



local circular velocity, Vc{r) — -y/— V<i> • r. That is, there 
are no random velocities, although the particles could 
be considered the guiding centers of a set of particles. 
Only one sub-halo is sent towards the ring in each sim- 
ulat ion. An individti al sub-halo is adequately modeled 
as a iHernquistI ()199Gf ) sphere with the mass-scale radius 
relation found in the Aquarius simulations, since for our 
purposes the exact details of the internal mass distribu- 
tion are not important. The sub-halo motion is simply 
at a fixed velocity in a straight line, so it crosses the ring 
only once which avoids the complications of a sub-halo 
orbiting and having multiple ring crossings. A single sub- 
halo is started at a distance of the ring radius from the 
ring. The time step is 0.001 unit of r/vc{r) at 30 kpc, or 
about 1.4 X 10^ yr. The ring of particles responds to the 
sub-halo and background halo potential but has no self- 
gravity. A series of simulations is done varying the angle 
with respect to the tangent to the ring in the direction of 
motion over = [41, 60, 90, 104, 120, 139] degrees, and 
the angle with respect to the plane of the ring oi 4> — [0, 
30, 60, 90, 270, 300, 330] degrees. Together these two dis- 
tributions non-redundantly coarsely sample an isotropic 
velocity dispersion field. The encounter velocities are 
varied over [0.3, 0.5, 0.8, 1.2, 1.6] times the circular ve- 
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locity. For an isotropic velocity dispersion of sub-halos in 
our potential the typical ring encounter velocity is close 
to Vc{r). The simulations are done for 10^, 10^, 10* and 
10^ Mq sub-halos. In total 9310 encounters are simulated 
for a ring. 

The outcome of a typical set of encounters after 6.92 
Gyr is shown in Figure [T] For purposes of easily vis- 
ible illustration we chose a relatively massive sub-halo, 
10* Mq, which has a scale radius of 0.83 kpc. The sub- 
halo is moving at 1.2vc. The top two rows show impact 
parameters of 0.06 kpc, where the sub-halo significantly 
penetrates the ring. The bottom row shows an encounter 
with an impact parameter of 1.5 kpc, about twice the 
sub-halo scale radius. The outcome at late times in all 
cases is that the ring is primarily perturbed perpendic- 
ular to the plane with loops pulled out. Note that in 
Figure [1] the vertical extent of the loops is stretched a 
factor of 220. The ratio of the upward and downward 
loop heights depends on the angle 0, as illustrated in the 
two upper rows. The upward or downward extension also 
varies with time as the perturbed particles pass through 
the orbital plane. For the large b relative to Rg passage 
of in the bottom row of Figure[l] the vertical deflection of 
the ring particles is somewhat larger than in the upper 
two rows, but there is essentially no visible gap in the 
ring density. The top two rows suggest that the density 
gap depth and length is insensitive to the angle (f>, the 
angle relative to the orbital plane. 

The gaps in these completely cold rings have a 
very simple morphology, whereas the gaps created in 
more realistic stellar streams with internal velocity dis- 
persion and a finite width are much more complex 
(|Carlberg et"aII[20TTI : lYoon. JohnstonlT Hogg 2011). In 
particular real gaps are often slanted to the centerline. 
Old streams can be subject to so many encounters that 
the remaining segments of the stream are relatively dense 
clumps smaller than the gaps between. If the progenitor 
is on a fairly epicyclic orbit the tidal stripping episodes 
occur at perigalacticon and lead to additional structure 
in the stream. Therefore the simplified model presented 
here is intended to be a statistical guide to the analysis 
of streams, not a detailed description of any particular 
stream. 

We measure the xy projected linear density relative 
to the unperturbed density using the distance, Adxy be- 
tween every third particle. The relative density is then 
nAdxy/(6T^r), where n is the number of particles in the 
ring and r is its radius. Figure [2] shows that the density 
profile, of the gap is effectively independent of </), the en- 
counter angle with respect to the plane. The particles 
that are stretched out in the loop pile up at the edges 
of the gap. The dependence of the depth of the den- 
sity dip in the gap with b is shown in Figure [31 showing 
that although a loop is pulled out, there is little density 
dip for encounters significantly beyond the sub-halo scale 
radius. 

2.1.1. Critical Impact Parameter 

The relative density at the bottom of the dip is dis- 
played as a function of the ratio of the impact parameter 
to the scale radius, b/Rs{M), in Figure HI The density 
measurements are averaged over all angles, 9 and (p. The 
results are displayed for velocities of 0.8, 1.2 and 1.6 Vc- 
Since the ring particles are moving at the circular ve- 
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Fig. 2. — The density distribution in the ring which a lO'' Mq 
sub-halo passed at right angles in the horizontal plane at a velocity 
of 1.2vc at fe = 0.3 kpc for the seven different angles, </> = [0°, 30°, 
60°, 90°, 270°, 300°, 330°], with respect to the orbital plane. 




O ^ — ' — I — ' — ' — ' — ' — I — ^ — ' — ' — I — ' — ' — ' — ' — I 

0.86 0.87 0.88 

fraction of circle 

Fig. 3. — The density distribution in the ring which a 10^ Mq 
sub-halo passed at right angles in the horizontal plane at a velocity 
of 1.2 times the circular velocity with a closest separation distance 
of 0.06, 0.6 and 3 kpc, which lead to increasingly shallow density 
minima. 

locity there is a range of encounter velocities with the 
particles in the ring. If the sub-halo velocity is signifi- 
cantly below the circular velocity there will be no very 
low velocity encounters in the frame of an individual star 
and the sub-halo, so on the average the gaps will be less 
deep. However since encounters below O.Svc comprise 
only about 10% of the encounters we leave them out of 
the average. 

We need to determine the maximum impact parameter 
that creates a visible gap. We will model the behavior of 
the density dip as a function of b/Rg (M) as a step tunc- 
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Fig. 4. — The minimum density of particles in the plane of the 
ring as a function of the distance of closest approach normalized to 
the scale radius of the sub-halo. The lines are for sub-halo masses 
of lO'*, 10'^, 10** and 10^ Mq from top to bottom, respectively. The 
measurements are made at 6.6 Gyr after the encounters for v = 
1.2vc. These measurements are made for 30 kpc radius rings. 



tion at a value a{M,r,t), with visible gaps created for 
b < aRs{M) and no gaps for larger h values. We evaluate 
the value of a = b/Rs{M) at | of the unperturbed den- 
sity. The resulting critical values of b/Rs{M) are func- 
tion of both mass and time, as shown in Figure [5l Using 
a density dip of i yields larger values of a, but its use 
turns out to give values within our overall modeling un- 
certainties. The function a{M, r, t) has a significant time 
behavior. At early times, roughly before 2 Gyr, the gaps 
are not well developed and we exclude them from the fit. 
Over a Hubble time the value of a nearly doubles. The 
time behavior is mainly relevant for dynamically young 
streams or the young part of streams that are still being 
created where we will expect the stream to be far less 
gapped than our average value. 
Fitting to the data of Figures |3] and O we find, 

a33%(M,r,t) - (1.71 + [0.28 + 0.101ogioM8][i-7]) 



(aokpc 



(2) 



where t is in Gyr and normally we will simply use the 
value at 7 Gyr as a reasonable average value for a{M, r). 
Fits which exclude the mildly non-linear 10^ Mq sub- 
halos give nearly equal values. 

2.1.2. The Gap Creation Rate 

We now have the elements to evaluate the gap cre- 
ation rate, Equation[T] The maximum impact parameter 
which creates a gap is bmax = a(M, r, t)Rs{M) where we 
will use the 7 Gyr value of a, roughly the average over a 
Hubble time. The mean approach speed to the ring is the 
integral over v±^ which is ^/Ij-na ~ O.SOcr for a Gaussian 
velocity distribution. The rate integral then becomes an 



Fig. 5.— The values of a(Af, r, i) = hjRs, at 1/3 of the unper- 
turbed density, at 30 and 60 kpc, solid and dotted lines respec- 
tively, for sub- halo masses of 10 , 10^, 10* and 10^ Mq from top to 
bottom, respectively. 

integral entirely in the mass distribution, 

f-OQ 

7^u(M,r) = 0.807r(Tn(r) / N{M)a33%{M,r)Rs{M)dM. 

(3) 

The Aquarius simulations (jSpringel et al.l[2008l ) find that 
Rs ~ 0.83M|-43 j^p(, a^j^d ^^^^^ ^he mass distribution func- 
tion is N{M)dM = 2MM^^-^dAh. The mass function 
gives the total number of sub-halos contained within 
the volume, V50, that has a mean interior density 50 
times the critical de nsity which has a radius 433 kpc 
(jSpringel et al.l [20081 ) . Putting these pieces together we 
find. 



7^u(A^,r) = 490 fiV; 



50 



"0 JMg 



M8"'-'a33(M8)Af8"-^^ dNh 
(4) 



where Afg = M/10 Mq is the smallest mass halo that 
creates a visible gap. The normalizing density in the 
V50 of the radial density variation, n(r), of the sub-halo 
population volume is no- The mass integral is dominated 
by the smaller masses so we have set the upper mass 
limit to infinity. Evaluating the integral, choosing a local 
a = 120 km s^^ and converting to astronomical units, we 
find. 



7^u(M,^■) = 0.0066 



lOOkpc 



n(r)/no 



6 120 kms"^ 
-°-35kpc-iGyr-i, (5) 

where we have used the iSpringel et al.l ()2008f ) measure- 
ment for n(r) (their Fig. 11) which finds that the sub- 
halo density at 100 kpc is 6 times the mean within the 
V50 volume. The rate at which gaps are created increases 
as the minimum mass sub-halo that can create a notice- 
able gap goes down, i.e. 7^ m (x M~^-^^ Ha lo-to-halo 
sub-halo number differences ()Ghen et al.l 120 111) will lead 
to ±50% or so variation in the rate in a specific halo. 
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2.2. Gap Loop Heights 

An important observable quantity is the length of the 
gaps. We measure the length as the distance in the or- 
bital plane between the two outside points at ^ of the 
vertical maximum of a loop. The length of a gap should 
be related to the distance that the loop of ring particles 
is pulled out of the orbital plane, that is, as the verti- 
cal extension of the ring increases, so does the length of 
the visible gap. This means we first need to understand 
the dynamics of the vertical perturbations. The impact 
approximation and the tidal approximation are the two 
general approaches used to describe the interaction of a 
satellite, such as a sub-halo, with a set of particles, the 
stars in the stream. The impact approximation is best 
suited to calculating deflections of individual stars from 
the rapid passage of a satellite and the tidal approxima- 
tion is best suited to problems involving the differential 
gravitational field across a finite sized object. Both ap- 
proaches are useful for calculating the gravitational heat- 
ing of the perturbed objects and at their root both have 
the same underlying gravitational dynamics. 

The impact approximation assumes that both a sub- 
halo and a star stream orbit can be approximated as 
straight lines and sets aside the possibility that orbital 
resonances play a significant role. A sub-halo of mass M 
moving past a star at a relative velocity v with a distance 
of closest approach b will induce a velocity perpendicu- 
lar to its direction of motion in the direction of closest 
approach of 6v± = 2GM/bv. Since many of the encoun- 
ters partially intersect the stream the relevant mass will 
be M{h), the mass enclosed within the impact parame- 
ter, not the total mass of a sub-halo. For an orbit that 
is approximately circular at radius r and taking 5v±^ to 
be in the vertical direction, then conservation of angular 
momentum requires that the maximum orbital deviation 
out of the plane is z„j — rSv± / Vc ■ 

We need to average over sub-halos crossing at all an- 
gles to the ring tangent and to the orbital plane. To 
a first approximation forces in the horizontal plane do 
not add to the angular momentum around the normal to 
the plane of the orbit, but do act to coherently increase 
the epicyclic motion around the guiding center, but that 
makes little difference to the stream density. A sub-halo 
moving at an angle to the orbital plane initially pulls the 
stream down and then up, which evolves into a paired 
set of peaks perpendicular to the orbital plane, as shown 
in Figure [T] The vertical extent of the perturbation is 
measured as the vertical distance from the positive peak 
to the negative peak, rather than the measurement of the 
absolute value of the offset from the initial orbital plane 
for a horizontal sub-halo passage. The resulting Z+ — Z~ 
distribution is shown as a binned average in Figure [6l 
The mean is about 30% higher than the prediction, with 
about a factor of two scatter above and below the line. 
For our use the scale offset is not an issue. 

The tidal approximation is an alternative description 
of the encounter of a sub-hal o with a stream. The Hill 
radius, Rhui — r-\/M/{3Mt), inside of which the gravi- 
tational field of the sub-halo, mass M, begins to dominate 
over the gravitational field of the total mass interior to 
the orbit of the stream, Mt . The plot of the vertical de- 
flections against Hill radius is shown in Figure [7] using 
the same data and same vertical scale as Figured Al- 
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Fig. 6. — The vertical extent at 33% of maximum height of the 
sub-halo induced "loop" to the impact parameter prediction for a 
ring at 60 kpc. The vertical extent and radial co-ordinate r are 
measured in scaled units of 30 kpc and the velocity in units of Vc. 
A ring at 30 kpc has half the vertical extent, as expected. 
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Fig. 7. — The maximum vertical height of a perturbation as a 
function of the tidal Hill Radius. The averages are slightly offset 
for clarity. 

though the vertical perturbation and the Hill radius are 
strongly correlated there remains a larger scatter than 
with the impact approximation. Including a power of b 
finds that RHiiib~^^^ removes some of the scatter. This 
is hardly surprising since with the addition of v this be- 
comes the impact approximation formula, multiplied by 
a constant. An alternative tidal quantity is the differen- 
tial tidal field, which is proportional to M/b^. A plot of 
the vertical deflections against M/b^ is mainly scatter. 
We conclude that the impact approximation is better at 
predicting the vertical deflections than the tidal approx- 
imation. 
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Fig. 8. — The ratio of the hnear size at 33% of maximum height 
of the sub-halo induced "loop" to the vertical height as a function 
of the impact parameter in ratio to the scale radius of the sub-halo. 
The lines show averages at fixed mass, starting with lO'' at the 
top increasing to 10^ Mq at the bottom. These values are for a 
ring at 60 kpc. 



2.2.1. The Characteristic Velocities 

Star streams have a mean particle motion close to the 
circular velocity. Therefore the encounter speed between 
a sub-halo and a star in the stream, which is required for 
the impact approximation, will on the average be higher 
than the mean circular velocity. A consequence is that 
low velocity encounters of sub-halos with the stream are 
extremely improbable. Taking the mean encounter along 
the stream to be Vc and the velocity dispersion in the two 
axes perpendicular to the stream as a the typical rela- 
tive velocity is ^Jvl + 2cr^. Using a ~ ■yc/1.8 in the outer 
part of a galactic halo then the typical sub-halo star rel- 
ative encounter velocity is \.21vc- The full triple velocity 
integral has the complication that there is a singularity 
at relative velocity of zero. Physically this would corre- 
spond to the velocity of the stream's progenitor sub-halo. 
Inserting a very small "core" velocity of about one me- 
ter per second, we find that the mean encounter velocity 
is l.QlTvc- Combining these factors we find that on the 
average, 

z (6) 

2.2.2. Gap Lengths 

The length of a sub-halo induced gap as a function 
of the sub-halo parameters is needed to determine what 
gaps are large enough to be visible. We measure the gap 
length at i of the maximum vertical height as a con- 
servative approach. The basis of a relation is visible in 
Figure [8] in which gap lengths relative to z,„ are plotted 
for a ring at 60 kpc. We see that for substantially pen- 
etrating encounters, b/Rs < 2, that the ratio of the gap 
length to the vertical extent, z^, = — , is roughly a 
constant at a given mass. Rings at smaller radii have rel- 
atively larger values of the length to height ratios, partly 



because Zm includes an r factor the ratio is dependent 
on the radius. Defining f3{M) as the length of the gap 
relative to the value of Zm we can find an approximate 
fit for the mass and radial dependence over the range of 
30 to 100 kpc to give. 



/333%(M,r) =60 



VSOkpc 



-0.63 



(7) 



The function /333% (M,r) appears to have no significant 
time dependence, once the gaps develop. At radii be- 
low 30 kpc the (3 value becomes somewhat larger than 
Equation [7] gives. 

The predicted length of a gap is the product of the 
average /3 and the impact Equation [5] for the vertical 
height of the loop in the star stream. 



Z(Af)=/333%(Af,r) 



2GM{Rs)r 
TmRjWji 



(8) 



Approximating a sub-halo density profile as a Hernquist 
profile we find th at the mass insi de the scale radius is 
M{Rs) = 1/4M HHernauistl [19901 ) ■ Evaluating the nu- 
merical factors we find, 



l{M,r) = 8.3 



30kpc 



0.37 



M°-4i kpc. 



(9) 



We now have an equation that gives the typical gap 
length as a function of sub-halo mass averaged over the 
velocity distribution. The gap length-mass relation can 
either be used as a consistency test, or, taken a step 
further and be used to eliminate the unobservable mini- 
mum mass parameter from the rate equation in favor of 
a minimum gap length. 

2.3. The Gap Rate Stream Width Relation 

The width of the stream and the minimum visible gap 
size are linked through orbital dynamics. The stream 
width is largely the result of stars emerging through the 
tidal Lagrangian points with low initial velocity disper- 
sion. The epicyclic approximation describes the subse- 
quent motion of the stars with small velocity dispersions 
around a guiding center. An epicycle has comparable 
amplitudes in the direction of the motion and perpen- 
dicular to the motion, hence the stream width and the 
random motion along the stream are of nearly equal size. 
Provided that the phase angles of the motion are not sig- 
nificantly correlated, any induced gap narrower than the 
stream width will be blurred out. Therefore we make 
the assumption that the width of the stream is equal the 
minimum gap size and use it to estimate the minimum 
mass sub-halo that creates a gap. 



w = 1{M). 



(10) 



Solving Equation [To] for M we substitute it into Equa- 
tion [5| to give the gap creation rate as a function of the 
width of stream, 



7^u(w,^) = 0.059 



n{r)/no 
6 



100 kpc 



0.55 



,-0.85 



kpc^^ 



for w in kpc, where we have used Equations [2] and [71 
Equation [TT] recasts the gap rate in terms of an observ- 
able, assuming that the CDM halo substructure model 
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Fig. 9. — The linear density of stars kpc"^ in the M31 NW star 
stream with moving average of 5.5 kpc (top) and 2.5 kpc (bottom). 
The hatched areas are the ±lcr average error of the mean. 

gives the correct density of sub-halos. Over the range of 
20-100 kpc where streams are generally fomid the factors 
of n{r) and r partially counteract each other to reduce 
the radial variation. Equation [TT] is an observational test 
of the CDM sub-halo prediction. 

3. GAPS IN OBSERVED STREAMS 

Counting gaps in a star stream requires considerable 
statistical care because the available observational data 
has streams embedded in a background which usually 
overwhelms the stream and is generally variable along 
the stream's length. A reliable measurement depends 
on having a good estimate of the density errors along 
the stream. The pseudo-images of stream star den- 
sities on the sky after filtering a star catalog with a 
color-magnitude relation can be usefully indicative but 
gaps should only be counted on extracted stream densi- 
ties with local error estimates. We concentrate on four 
streams, the NW stream of M31, and the Pal 5, East- 
ern Banded Structure and Orphan streams of the Milky 
Way. All of these have published estimates of the linear 
density of stars along the stream. 

3.1. The North- West Stream of M31 

The overall M31 PAnd AS imaging project a nd th e 
streams ar e described in iMcConnachie et"al] ()2009D : 
iRichardson et al. (201l|) and t he details of the NW 
stream in TCarlberg et all (poTl . Here we use the lin- 
ear density profile to estimate the number of gaps. In 
Figure IHl we show the density field as a function of the 
angle around the s tar stream, with all q uantities mea- 
sured as reported in lCarlberg et al.l ()2011[ ). The plots are 
made with a moving average over 5 and 11 bins of 0.25°, 
or roughly 2.5 and 5.5 kpc, respectively. To determine 
what scales in the stream appear to be physically corre- 
lated we plot the cross-correlation of the stream density 
minus the mean value with itself in Figure [TUl We see 
that the stream densities are positively correlated over 
the first three bins of 0.25°, after which the fiuctuations 
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Fig. 10. — The cross-correlation of the density field measured in 
0.25° bins. The correlation is positive for the first 4 or 5 points, 
indicating a scale of about 1 degree. 



are essentially random. This correlation length is com- 
parable to the width of the stream and helps validate 
the assumption that leads to Equation[TUl Therefore the 
smallest size scale we search for gaps in these data is 
about 1°, or, 2 kpc. 

To count the gaps we need to have a good measure of 
the level of statistical fiuctuations. The cross-correlation 
of the data with itself, at shifts beyond 5 bins, where 
there is no intrinsic signal, provides an accurate inter- 
nal estimate of the variance of the density field. The 
product of two normally distributed random variables is 
distributed like a modified Bessel function of the second 
kind and has a large kurtosis, 9, as compared to a nor- 
mal distribution with a kurtosis of 3. However, the sum 
of many products (115 in our 0.25° sample) converges to 
a Gaussian distribution in accord with the central limit 
theorem. The kurtosis of n sums of products is 3 -I- 6/n 
so declines fairly quickly to the Gaussian value of 3 with 
increasing n. The variance around the mean density of 
the data within the stream region, prior to background 
subtraction, is 1.22 (stars/kpc)^, and a comparable 1.33 
in the background region. The background subtracted 
mean stream density has a variance of 6.16 (stars/kpc)^. 
The F ratio of the measured variance over the randomly 
expected variance is 6. 16/(1. 22-1-1. 33)=2. 42, which has 
a random chance of occurrence of less than 2 x 10~^. 
Larger bins give even higher significance that the mea- 
sured stream density has variations well in excess of what 
is expected from random fluctuations. The conclusion 
and the levels of s ignificance are i n acco rd with the re- 
sults presented in iCarlberg et al.l ()2011[ ) but use an in- 
dependent line of argument which does not rely on an 
explicit calculation of shot noise statistics. 

An estimate of the number of gaps is simply to count 
the number of times the density crosses down and up 
through some density threshold, requiring that the gap 
be at least the size of the smooth length. For a 5 point 
moving average at [— Icr, 0, -|-1(t] density levels we find 
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[11,11,15] gaps of mean length [7.1,6.9,15.1] bins, or, [1.8, 
1.7, 3.8] kpc. Taking the stream length to be 200 kpc 
and assuming an age of 10 Gyr we find from the 12 ± 
2 gaps that the M31 NW stream has been subject to 
a mean gap creation rate of 0.006 ± 50% gaps kpc~^ 
Gyr~^. In principle the predicted rate also includes any 
sufficiently massive object, but the detectable mass is 
well above the range for globular clusters and we do not 
expect molecular clouds along the large radius orbit of 
the NW stream. 

The rate equation. Equation [SJ applied to the M31 
NW stream requires a minimum sub-halo mass of M = 
1.3 X 10* M0 and implies a total of 205 sub-halos, which 
is a fac tor of 7.3 higher than 28 known dwarf galaxies 
of M31 (i Richardson et ai1l2011h but includes objects out 
to 433 kpc, of which the only the inner 150 kpc radius 
volume has been carefully searched. Therefore the M31 
stream not demand a huge dark matter sub-halo pop- 
ulation by itself, but this analysis does recover the well 
known discrepancy which begins to appear for even fairly 
massive sub-halos. 

3.2. The Pal 5 Stream 

A narrow stream is predicted to have more visible gaps 
relative to wider stream at larger radius. The globu- 
lar cluster Pal 5 has both leading and trailing streams 
with the northern, trailing, arm having the most data 
in the Sl oan survey. With data made available be- 
fore DRl lOdenkirchen et all ()2003l ) found the stream 
to be about 6 degrees long and measured the width to 
be 0. 11 kpc. With the better c overage of the SDSS 
DR4. lGrillmair fc DionatosI (|2006[) found that the north- 
ern arm was about 18.5° and provided a surface den- 
sity contour plot showing the considerable variation in 
the surface density of the stream. The density varia- 
tions closest to the cluster are convincingly explained as 
epi cyclic pileups of stars tidally stripped from the clus- 
ter iKiipper et al.l (|2010l 1201 ID . however beyond about 
Ikpc (2.5 degrees) the variation in phase angles smooths 
out this behavior. F rom the published Figure 3 of 
iGriUmair fc DionatosI (poM ) we count regions in the 
stream below the lowest contour to find that there are 
about 5 gap s over the ou ter 16° of the northern tail of 
Pal 5. lOdenkirchen et al.l (20031) estimated that the 5-6° 
of stream took about 2 Gyr to fill out, which scaled to 
the 18.5° length gives an age of about 7 Gyr. With these 
values we estimate the linear gap creation rate in the 
northern stream of Pal 5 to be about 0.11 kpc~^ Gyr~^ 
to which we assign an error of 50% based on uncertainties 
in the gap counts and stream age. 

3.3. The EES Stream 

iGrillmairl ()2011| ) finds that the Eastern Banded Struc- 
ture is 0.17 kpc wide and located at a galactocentric dis- 
tance of about 15 kpc. Figure 4 of his paper plots the 
density along the 18° of the stream, with an estimate 
error of 30% in each bin of 0.5°. We estimate that there 
are between about 4 and 12 gaps along the stream. We 
will take the stream age to be approximately 7 Gyr to 
derive a gap creation rate of 0.24 gaps kpc~^ Gyr~^, with 
about a factor of two uncertainty in this value. 

The estimated rate of gap creation in the EBS stream 
with Equation [5] implies a minimum halo mass of M ~ 



TABLE 1 
Observed Stream Gap Statistics 



Stream 


Gaps 


Length 


Width 


Age 


Hoc 


n/no 




# 


[kpc] 


[kpc] 


[Gyr] 


[kpc] 




M31-NW 


12 


200 


5 


10 


100 


6 


Pal 5 


5 


6.5 


0.11 


7 


19 


24 


EBS 


8 


4.7 


0.17 


7 


15 


30 


Orphan 


2 


30 


1.0 


3.9 


30 


20 



1 X IO^IvIq. The derived minimum mass implies a total 
of 1 X 10^ sub-halos, although we note that both the 
mass and in turn the implied total number of sub-halos 
are sensitive to the gap rate estimate and the details of 
the model in Equation [5l Possible complications for the 
EBS stream is that the sub-halos are subject t o erosion 
by the disk and bulge ()D'Onghia et al.l 120101 ) and the 
stream may encounter smaller scale structures such as 
spiral arms and molecular clouds which could add to gaps 
in the stream. 

3.4. The Orphan Stream 
Th e Orph an Stream (|Grillmairi 120061: IBelokurov et al.l 



120061 12007D has received considerable attent ion, culmi 
nating with a detailed orbital model ((Newberg et al.l 
l2010f ). In creasing the me an distance from 20 to 30 
kpc revises iGrillmaiii (|2006D 's estimated width of 700 pc 
to 1 kpc. Gaps are marginally present in Figure 5 of 
iNewberg et al.l (|2010 ') at -16° and -31°, particularly when 
compared to the bottom panel of their Figure 17. We es- 
timate that there are 2 gaps prese nt, although we will al - 
low a 100% error on this estimate. INewberg et"all (|2010f ) 
present a model that estimates the age of the stream to 
be about 3.9 Gyr. If we had used the dynamical age 
of 3.9 Gyr rather than a uniform 7 Gyr the correction 
would raise the Orphan data into better agreement with 
the prediction. 

4. GOIvIPARISON OF PREDICTIONS AND 
OBSERVATIONS 

The observed stream gap rate and width estimates are 
summarized in Table 1. The results of this paper are 
largely encapsulated in Figure 1111 which plots the mea- 
sured and predicted gap rates against the widths of the 
streams. The data have all been rescaled to a galac- 
tocentric distance of 100 kpc using n{r)/no and r in 
Equation [11] The agreement of the data and the CDIVI 
sub-halo prediction is remarkable. We show the relation 
for an arbitrary shallower N{M) oc M~^-^ normalized 
to have the same number, 33, sub-halos above 10^ IVIq 
present in the LCDM prediction. If one used the vis- 
ible dwarf galaxies to predict the number of sub-halos 
the predicted density of gaps in these streams would be 
roughly a factor of 100 smaller, which would be in wild 
disagreement with the observation of virtually any gaps 
at all. The numerical consistency that we find is both an 
indirect and statistical test for the existence of the large 
predicted population of dark matter sub-halos. 

The Orphan stream merits further observational at- 
tention since the available measurements are potentially 
consistent with it being a stream with no gaps. There 
is room for variation in these statistical estimates, but 
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width [kpc] 



Fig. 11. — The estimated gap rate vs stream width relation for 
M31 NW, Pal 5, the EBS and the CDM halo prediction. All data 
have been normalized to 100 kpc. The width of the theoretical re- 
lation is evaluated from the dispersion in the length-height relation 
of Fig. \S\ Predictions for an arbitrary alternative mass functions, 
N{M) Qc Af~^'® normalized to have 33 halos above 10^ Mq is 
shown with a dotted line. 



when the other streams work fairly well with the adopted 
CDM spectrum the Orphan stream is expected to fall 
into line and could potentially be a case against a rich 
sub-halo population in our galaxy. 



5. DISCUSSION AND CONCLUSIONS 

This paper lays a basis for understanding the rate at 
which dark matter sub-halos create visible gaps in stellar 
streams. Restricted n-body experiments guide a calcula- 
tion based on the impact approximation. The outcome 
is a prediction of the local rate at which gaps are in- 
duced in a stellar stream which effectively counts the 
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number of sub-halos massive enough to produce visible 
gaps. The rate of gap creation is cast as a relationship 
with the stream width, which gives a predicted relation- 
ship between observables that tests the CDM sub-halo 
prediction. 

Comparison of the CDM based prediction of the gap 
rate- width relation with published data for four streams 
shows generally good agreement within the fairly large 
measurement errors. The result is a statistical argument 
that the vast predicted population of sub-halos is indeed 
present in the halos of galaxies like M31 and the Milky 
Way. The data do tend to be somewhat below the pre- 
diction at most points. This could be the result of many 
factors, such as the total population of sub-halos is ex- 
pected to vary significantly from galaxy to galaxy, allow- 
ing for the stream age would lower the predicted number 
of gaps for the Orphan stream and possibly others as 
well, and most importantly these are idealized stream 
models. 

There are many improvements possible to strengthen, 
or, destroy this result. Each observed stream could be 
realized through an n-body model and inserted into an 
appropriate model of a sub-halo rich galactic halo. The 
result would still be statistical but the various time de- 
pendent and mass dependencies would be fully taken into 
account. On the observational side the presence or ab- 
sence of a particular gap in a stream are somewhat uncer- 
tain due to the relatively low signal to noise of streams 
at the present time. A statistically robust measure of 
the number of gaps needs to be developed. However, 
the observations are improving quickly and much more 
secure gap descriptions will soon be available. Overall, 
within our simplified general modeling of sub-halos cross- 
ing streams to create gaps, we finding a perhaps surpris- 
ingly good agreement between data acquired for other 
purposes and the CDM prediction of a large population 
of sub-halos. 
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